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Abstract 

The conspicuous threshold enhancement has been observed in the baryon-antibaryon subchannels 
of many three-body B decay modes. By examining the partial waves of baryon-antibaryon, we first 
show for B^ — > ppK that the pK^ angular correlation rules out dominance of a single pp partial 
wave for the pp enhancement, for instance, the resonance hypothesis or the strong final-state 
interaction in a single channel. The measured pK angular correlation turns out to be opposite to 
the theoretical expectation of a simple short-distance picture. We study the origin of this reversed 
angular correlation in the context of the pp partial waves and argue that NN bound states may 
be the cause of this sign reversal. Dependence of the angular correlation on the pp invariant mass 
is important to probe the underlying issue from the experimental side. 
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I. INTRODUCTION 



In the baryonic B decay the three-body modes dominate over the two-body modes. Fur- 
thermore, in the three-body decay,the baryon-antibaryon pair is copiously produced at small 
invariant mass near the threshold [l|, H, 0, Q, [5[ . Various theoretical ideas p, Q, H, 0, 13], some 
kinematical and others dynamical, were proposed for this threshold enhancement. A simple 
short- distance (SD) argument can explain qualitatively both the dominance of three-body 
modes and the threshold enhancement of baryon-antibaryon: To produce a baryon and an 
antibaryon in the two-body decay (Fig. la), one energetic qq pair must be emitted back 
to back by a gluon so that the gluon emitting the qq pair is highly off mass shell. The 
hard off-shell gluon suppresses two-body decay amplitudes by the power of a s /t, where t 
is the four-momentum square transferred through the gluon. In the three-body decay with 
an additional meson (Fig. lb), a baryon-antibaryon pair can be emitted collinearly against 
the energetic boson in the final state. In this configuration a quark and an antiquark are 
emitted by a gluon nearly in the same direction so that the gluon is close to the mass shell 
and the short-distance suppression does not occur. Consequently the pp of small invariant 
mass is strongly favoured. 





FIG. 1: Short-distance picture in quarks and antiquarks for (a) two-body baryonic decay and (b) 
three-body baryonic decay. In the two-body decay (a) the fat virtual gluon (the thick vertical solid 
line) must split into qq while in the three-body decay (b) nearly on-shell gluons (broken lines) turn 
into qq. The slow spectator antiquark is denoted by the short line q s . 

In addition to threshold enhancement, the angular correlation was measured between the 
final proton or antiproton and the boson in some modes, most clearly in B ± — > ppK^ (ll| . 
Then an intriguing puzzle[Io[ has emerged in the preceding SD picture: In that picture, 
the antiproton momentum should point more likely to the direction of the K~ momentum 
in the pp rest frame of B~(bu) decay. That is, the proton should tend to move away from 
K~ in this frame. The reason is that the antiproton p picks up the slow spectator w-quark 
and therefore its momentum is smaller on average than that of the proton p in the B~ rest 
frame. By boosting the B~ rest frame to the pp rest frame, we reach this conclusion. 



However, the Belle Collaboration showed exactly the opposite[TJj; it is the proton that 
is emitted along K~ in the pp rest frame. Belle selected the threshold events by making a 
cut in the pp invariance mass m P p{< 2.85 GeV), but did not give the angular correlation as 
a function of m pp - for the selected events. Meanwhile BaBar gave a Dalitz plot of ppfT~[3| 
from which one can read the same trend as Belle's angular dependence. 
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Rosner|8j argued qualitatively in terms of quark diagrams 12j and predicted this angular 
correlation with baryon production through diquarks. But the ar gum ent does not seem to 
work for all baryonic modes in its simple form. Cheng and others jlOf] computed the d ecay 
amplitudes in the pole model with factorization, leaving out inelastic form factor terms 0, 10 . 
Their result does not lead to the correct angular correlation in the case of B~ — > ppK~ . 
The simple SD picture presented at the beginning is successful in the angular correlation of 
most three-body baryonic modes, e.g., B~ — > Ap7, but fails notably for B~ — > ppK~. Its 
failure suggests us importance of long-distance (LD) effects somewhere in the decay process. 
Indeed, the fragmentation by quark diagram and the pole model both contain some of the LD 
effects in very different ways. In this paper we take a close look at this angular correlation of 
B ± — ► ppK ± from the viewpoint of partial waves in general and try to resurrect the simple 
SD picture by incorporating an appropriate LD effect in it. 

In our proposed analysis we first examine the partial- wave content of pp in B~ — > ppK~ 
(and its conjugate) and conclude purely kinematically that the pp enhancement cannot be 
a broad resonance. For the same reason we rule out strong final-state interaction (FSI) in 
a single pp partial wave as a cause of the enhancement. We shall observe that reversal of 
the angular correlation occurs if some LD effect flips relative signs of partial-wave decay 
amplitudes. Such sign flip may indeed occur if NN bound states exist in right channels. 
The recently discovered state X(1835)[13] is a good candidate that may be responsible for 
the sign flip. If X(1835) should be an NN bound state, we expect a similar bound state in 
other channels from our reasoning of the sign flip. 



II. pp PARTIAL WAVES IN B~ -» ppK~ 

We study the angular correlation between the proton momentum and the kaon momentum 
in the rest frame of pp by choosing the z-axis along the K~ momentum. (Fig. 2) 




9s* K 



FIG. 2: The pK angular correlation in the pp rest frame. 



Since B~ meson and K~ meson are spinless and their momenta are both along the z- 
direction in the pp rest frame, the ^-component of total angular momentum is zero for pp 
by AJ Z = in this frame. Following the standard helicity formalism 141]. we can describe 
the angular dependence of the helicity decay amplitudes with Wigner's ci-functions [15] as 



A(B -> ppK~ 



E ^A p A f ;0^0a( 



-i\4>p 



(A 



(1) 
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where X p and Xp are the helicities of p and p in the pp rest frame, (9 P , 4> p ) are the angles of the 
proton momentum in this frame (Fig. 2), and A x ^ x _. is a function of the pp invariant mass 
mpp. Since experiment does not measure helicity of proton nor antiproton but sums over 
all helicity states in what follows, the differential decay rates are <j) p independent. Therefore 
we need not specify the direction of P = in our case. 2?J Squaring the amplitude and 



summing over the pp helicities X p and Xp, we obtain the differential decay rate: 



dT(6 p ) 



r, 



B -^ PP K- XpX _ 



o£|£<V&(*p) , (2) 



where T includes kinematical factors that depend on m P p. If we make the usual assumption 
that the strong penguin interaction dominates in the decay B — > ppK, the CP-violating 
phases drop out of the decay rate. Under parity reflection the angle 8 p remains unchanged, 
while under charge conjugation the angle 9 P turns into tt — 9 P of B + — ► ppK + because of the 
interchange p <->p and K~ «-> K + . Therefore, 



dT(9 p 



diTDjipipd^hip 



dT(n - 9 p ) 



B+^ppK+ arripp 



(3) 

B-^ppK- 



We shall be able to use this equality as a test of the penguin dominance. The corresponding 
relation holds between B° — > ppK and B — > ppK . 

At this stage we can prove that the pp enhancement is not a resonance, for instance, a 
glueballQ: The Wigner functions dQ X (9), which are proportional to the associated Legendre 
functions, possess a special symmetry property under 9 <-> 7r — 0[H|, 

d J 0X (n - 9) = (-l) J ^d J 0X (9) -> I^aC^-^I 2 = I^aWI 2 - (4) 

If the pair is produced entirely through a resonance, only the term of the resonance spin 
J contributes in Eq. (j2j) without sum over J. Since the function |dg A (#)| 2 is unchanged 
under 9 — > it — 9 (i.e., cos 9 — > — cos#), so is dT/dfl p in this case. However, experiment 
shows a pronounced asymmetry between two hemispheres of cos > and cos < 0. (Fig. 



3.) In terms of the forward-backward asymmetry parameter 11 



A=(N+- N-)/(N + + AL) = 0.59™ (5) 

contrary to A = in the case of a single J. Although interference between the resonant 
and nonresonant amplitudes of different J can produce some asymmetry in principle, such 
interference should be insignificant under the normal circumstance where the resonant ampli- 
tude acquires the phase | through the resonant decay relative to the nonresonant amplitude: 
arg(A^ es A^* n ) ~ ±7r/2. It is implicitly assumed here as usual that the nonresonant produc- 
tion amplitude dose not acquire a significant phase. If for some reason the large asymmetry 
of Fig. 3 should be caused by the interference between the resonant and nonresonant ampli- 
tudes of different J's, the very small yield observed toward 9 P = n in Fig. 3 would mean 
nearly perfect destructive interference between them. In this case the nonresonant yield 
would have to be just as large as the resonant one. Therefore the marked asymmetry in the 
angular correlation rules out convincingly the hypothesis of pure resonant production. If one 
attempts to explain the the enhancement by strong FSI in a single dominant partial-wave 



channel of pp[16|. one would likewise obtain A ~ for the angular correlation. To be consis- 



tent with the observed angular correlation, partial-wave amplitudes of even and odd J must 
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FIG. 3: The pK~ angular distribution in the pp rest frame [Ref. 10]. 

coexist and almost maximally interfere. Our argument is very general and independent of 
dynamics up to this point. We now proceed to take dynamics into account. 



In the SD picture the spectator w-quark of B~ enters the antiproton with two energetic 
antiquarks (u and d) which are pair-produced nearly collinearly by two gluons. (Fig. lb) 
Note that it is a color-suppressed process for the s-quark from b — ► sg* of the strong penguin 
decay to form K~ by capturing the spectator ?Z.{28j In contrast, the proton consists of three 
energetic quarks; one from the primary decay interaction and two of pair-produced quarks. 
In the B~ rest frame, therefore, the proton recoils against K~ more energetically on average 
than the antiproton does. It means that in the pp rest frame, the proton tends to move away 
from K~ faster than the antiproton does. That is, A < contrary to the measurement. 
This is the "angular correlation puzzle" . No reasonable explanation has been given from the 
SD viewpoint. We must look for some LD interaction effect that has not been commonly 
appreciated. 

The maximum of the pp enhancement occurs near m pp - = 2 GeV in the BaBar data [3] and 



roughly <2.2 GeV in the Belle data[ll(. The dominant relative orbital angular momenta 
of pp are expected to be s-wave and p-wave. The amount of ci-wave is presumably small 
and higher waves are even smaller. The terms that contribute dominantly in Eq. (pQ) are 
therefore J = (X^o), J = 1 ( 3 Si, 3 PiM), and J = 2 ( 3 P 2 )- The explicit forms of the 
relevant d$ x functions (J < 2, A = —1, 0, +1) are[l5|: 

d° 00 (6) = 1, <4(0)=cos0, d 2 00 (6) = (3 cos 2 6- l)/2, 



4±i(0) = TV 1 / 2 sin 0, ^ ±1 (^ = ^3/2sin^cos^ (6) 

It is convenient to rearrange the helicity decay amplitudes A J XpX _. Q with the spectroscopic 

notation into A( 2S+1 Lj] X p — Xp). When only s-waves and p-waves are retained, the helicity 
amplitudes of definite isospin / for iViV can be written as 

A° 7 ±;0 = ±A I ( 1 S ,0)+A I ( 3 P ,0), 
A" = A I ( 3 S U 0)±A I ( 1 P 1 ,0), 
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A 11 



V2A I { 3 S 1 ;±1)±A I { 3 P 1 ;±1) 
A 2 * = A'( 3 P 2 ;0), 



.4 
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±=F;0 



3/2A J ( 3 P 2 ; ±1) 



(7) 



where we have denoted the helicity indices A p , Xp = ±1/2 simply by ±. We shall use 
this notation hereafter. The normalization of the amplitudes is arbitrary for the decay 
amplitudes which have no unitarity constraint. Since A J = for the strong penguin decay, 
the decay amplitudes for the charge eigenstates of B and K are given by the decay amplitudes 
A I ( 2S + 1 L J , A) of definite NN isospin as 



.4 



K~ (2S+1 



A 



K (2S+1 



Lj,X) 
Lj,\) 



1/2 A 1 ( 2S+1 L J , A) - A°( 2S+1 Lj, A) 
1 < 2S+1 L J ,\) + A°e s+1 Lj,\) 



1/2 



.4 1 



(8) 



Combining Eqs. (j7j) and (jSj), we obtain the decay amplitudes as functions of 9 p , 4> p and m pp . 



-ftO 



The decay amplitudes for B + /B° — > ppK + /ppK° are obtained from those of B /B 



ppK~ jppK with the interchange A p <->• Xp followed by 6 P — > 7r — # p and P — > 7r + p up to 
the overall CP phase factor of the penguin decay. 

We are now able to write the complete differential decay rate for B~ — > ppK~ with in 
s and p- waves in the notation of 23+1 Lj\ 



dY 



(A R - ('So, 0) + A K ~ ( 3 P , 0)) + (A K ~ ( 3 S u 0) + A K ~ CP,, 0)) cos 6 p 



+ ^~( 3 P 2 ,0)(3cos 2 # p -l)/2 



+ 



(—A K ~ CS , 0) + A K ~ ( 3 P , 0)) + (A k ~ ( 3 S u 0) - A K ~ ( X P 1; 0)) cos # P 



+ ^~( 3 P 2 ,0)(3cos 2 fl p -l)/2 



+ 



(A K ~( 3 S 1 , 1) + Jl/2A K -{*P X , 1)) sin# p 



+ A 



P 2 ,l)(3sinfl p cosfl p )/2 



+ 



1/2A 



+ ( 3 P 2 ,-l)(3sin# p cosfl p )/2 



Pl, -1)J sin# p 

2 



(9) 



where the P dependence goes away from the squared amplitudes of definite helicity A. 
Before going further, we point out that the s-wave amplitudes alone cannot generate the 
asymmetric angular correlation for B~ —> ppK~ even though two s-wave amplitudes {}Sq 
and 3 S\) enter the right-hand side of Eq. (Q: The reason is that the interference terms 
cancel out between 1 Sq and 3 Si as 



dT 



d/TDjrp'pd^h'p 



B-^ppK- 



T [\A K ( L S ,0)+A K ( 3 S u 0)cos9 p 
| - A K ~ CSq, 0) + A K ~ ( 3 S!, 0) cos 9 P \ 



+ |A^( 3 5 , 1 ,l)| 2 sin 2 ^ + |A K_ ( 3 1 S' 1 ,-l)| 2 sin^ ; 



(10) 
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and consequently dT/dfl p turns out to be symmetric under cos9 p — > —cos9 p . The same 
statement holds valid for p-waves alone. The observed steep asymmetry (Fig. 3) requires 
more than one orbital angular momentum, most likely s-wave and p-wave. It is very im- 
portant experimentally to study how the angular correlation varies as p-waves increase with 
m pp relative to s- waves across the threshold enhancement. It does not make sense to make 
a theoretical fit to the shape of the m pp plot without large interference between different pp 
partial-waves. 

The experimental uncertainty in the angular measurement limits quantitative analysis 
at present. Let us be content with qualitative analysis in this paper by approximating or 
interpreting for simplicity the angular correlation in Fig. 3 as ~ (1 + cos6 p ) 2 . This cos6 p 
dependence is realized if 

^-(^o.O) ~A K -CS!,0) ~A K ~CPi,0) ^ A K -( 3 P ,0), (Exp) (11) 

and all other amplitudes are negligible. An alternative solution is 

A K -CS ,0) ~ -A K -( 3 S 1 ,0) ~ A K -( 1 P 1 ,0) ~ -A K -( 3 P ,0), (Exp) (12) 

and all others are negligible. A small amount of A K ( 3 P2, 0) with the same sign as 
A K ( 3 Si, 0) would improve the fit a little by lowering the curve near cos8 p = and raising 
it near cos8 p = ±1, but it is not crucial to the essence of our qualitative argument. While 
an accurate prediction is difficult because of our deficiency in knowledge of the quark dis- 
tribution in baryons, the SD argument predicts, as we have argued above, the sign of slope 
opposite to Fig. 3: The angular dependence should be more like (1 — \a\ cos9 p ) 2 (\a\ < 1) 
in the SD argument. This angular dependence corresponds to the partial- wave amplitudes, 

A K ~ C$0,0) « A^(XO) « -A K ~{ 1 P l ^) « -A K '-( 3 P ,0), (SD) (13) 

or alternatively, 

A K ~ C$0,0) « -A K ~( 3 S l ,0) « -A K ~CPi,0) « A R -( 3 P ,0), (SD) (14) 

instead of Eqs. or (|12[) . Comparing the SD prediction with experiment, we find that 
the relative signs of the s-to-p-wave amplitudes are opposite. There are several alternatives 
that can alter the SD prediction in line with experiment: Sign reversal of the 1 Sq and 3 S± 
amplitudes brings Eq. ({TBI to Eq. (fill and Eq. (Till) to Eq. ([12]) . Alternatively, sign 
reversal of 1 Sq and 3 Pq brings Eq. ({TBI to Eq. ([121 and Eq. (fill) to Eq. (TTTT) . Sign reversal 
of 3 £i and l Pi also accomplishes the same. We ask what LD effect can possibly cause the 
sign reversal from Eq. (I13p or (1141) to Eq. (1111) or (IT2|) . In the next section we argue that 
the desired sign reversal may occur if bound states exist in some of the pp channels. Unlike 
the argument that has ruled out a pp resonance, this is speculative and admittedly less clean 
part of our argument. 



III. FINAL-STATE INTERACTION 



The three-body final-state interaction FSI was analysed in the approximation of sum 



of two-body FSI since going beyond is mathematically formidable 17]. Fortunately, in the 



particle configuration of our interest where the invariant mass of pp is small and the K 
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meson recoils fast against pp, it is a good approximation and at least a common practice 
to separate the two-body FSI of pp ignoring the rest of FSI. Inclusion of pp annihilation 
channels is more a difficult problem. If one wants to make a quantitative analysis, this will 
be a main source of uncertainty. [2^] Our task here is not to obtain numerically accurate 
results but to search a possible cause of sign flip for the amplitudes in the FSI. In order to 
make the sign flip argument plausible, we do not need much more than a basic argument of 
the elastic two-body FSI and its diagrammatic explanation. 

The standard practice in FSI resorts to potential theory and incorporates FSI by modi- 
fying the decay amplitudes with final particle rescattering as[2l| 

(15) 



A JI {s) — A JI (s)/f JI (-k) 



Jii 



where f JI (k) stands for the Jost function |22| of a partial-wave eigenchannel in variable 

k = \ J~s — Am 2 N (s = (p p + pp) 2 ). It is normalized to f JI (oo) = 1. This FSI factor sums 
up ladders or bubbles of final particle rescattering in potential. The Jost function can be 
expressed with the phase of scattering amplitude 5 JI 



in the Omnes representation 23 



f JI (~ 
A JI 



v_ 

7T 



6 JI (i/) 



"0 



v' — V 
' S JI ( 



V 



-dv', 

-du' + i5 JI ( 



v — xe 



(16) 



where v = k 2 
region (s < 4m 2 A 



The lower bound u of the dispersion integral is extended to the negative 
N j when pp annihilation into meson channels is taken into account. 



If annihilation and inelastic scattering are ignored, the phase 5 JI (u) would be equal to 
the phase shift of NN scattering according to the so-called Watson's theorem 24j . If there 



is a resonance in this elastic case, the phase shift 8 (u) rises from zero, passes through 7r/2 
at the resonance [y = vr) and approaches 7r as v — ► 00. (Fig. 4.) Therefore the phase of 
the decay amplitude acquires a minus sign (= e l7T ) above the resonance v = vr. 



5(v) 




FIG. 4: The phase S(u) of the FSI factor e A// across an elastic resonance. 



This negative sign is easily understood in diagram. When a final particle pair is produced 



through a resonance, as depicted in Fig. 5 for B 
resonance takes the form of 



ppK , the decay amplitude near the 



.4 



Jl 1 



A 



j 1 , 



9 2 (s) 



m\ - im R T R (s) - s' 



(17) 
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where A JI (s) is the amplitude in the absence of a resonance, m R and Y R {m 2 R ) are the 
resonance mass and width, and g 2 {s) is positive at s — m 2 R . The reason for positivity of 



g 2 {m 2 R ) is as follows: By the phase theorem of FSI 24j, it holds that &rg[A JI (s) / A JI (s)] = 
S JI (s). Therefore, the phase of A JI (s) / A U \s) must be equal to +vr/2, not —n/2, at the 
resonance peak. The phase of the resonant FSI factor in Eq. (|17p must agree with this value 
at the resonance. The phase of aig[ig 2 j 'm^r(m^)] is +7r/2 at s = m|, in agreement with 
the phase theorem if g 2 {m 2 R ) > 0. The function g 2 (s) is expected to be only mildly energy 
dependent even near the threshold since the centrifugal factor k l of orbital motion resides 
in A JI (s) and A (s), not in g 2 (s). Since g 2 (s) (~ g 2 {m R )) is positive, the sign of A JI (s) 

is the same as that of A J \s) below the resonance (m 2 R — s ^> m#r#), but turns opposite 
above the resonance (s — m 2 R 3> m R Y R ). This is another way of seeing the sign and energy 
dependence of the FSI factor across a resonance. 




FIG. 5: Resonant production of pp in B — > ppK . 

This simple argument is modified by inelasticity above the NNir threshold and by anni- 
hilation into meson channels. Above the energies where inelastic channels start contributing 
substantially at ^Js > 2m at + m n , the FSI formulas of potential theory is no longer ap- 
plicable. If we simply truncate the phase integral in Eq. (I16p at v — z/ max somewhere 
above the inelastic threshold, the FSI factor computed in the narrow-width (step-function) 
approximation turns out to be 

A JI (u) ^max ~ ^ / _^ \ /10 n 

e , (z/>z/ R ), (18) 

v R - v 

which satisfies f(—k) ^ 1 as v[= k 2 ) — > oo. This FSI factor is negative between the 
resonance and the inelastic threshold; 

l/f(-k) = e AJ1 ^ < 0, (y R <v< u max ). (19) 

It means that the FSI factor gives a minus sign above the resonance until energy goes up 
so high that inelasticity becomes important. While the negative sign is easy to understand, 
magnitude of the FSI factor is harder to estimate since it depends on the dispersion integral 
over the entire energy range. 

Let us turn to the effect of the annihilation channels into mesons. The first issue is that 
the phase 5 JI of the decay amplitude is no longer equal to the phase of the NN scattering 
amplitude at any energy where annihilation occurs. An approximate equality between two 
phases holds only in those eigenchannels in which the channel coupling is weak between iViV 
and the meson channels. It is not obvious whether this is the case for the relevant iViV 
channels near the threshold. We must assume it here. Discussion will be made on this point 
below and in the next section. The other issue is whether relevant NN resonances really 
exist or not. The candidates of NN bound states and resonances indeed exist. Since the 
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NN bound states can be only loosely bound, N and N are spatially separated outside the 
range of annihilation interaction and therefore the annihilation into mesons is suppressed. 
Meanwhile, being a bound state, the state cannot decay into a nucleon and an antinucleon 
since its mass mg is below the NN threshold (tub < 2m^)- However, there is an escape 
from this argument: The finite lifetime due to meson annihilation generates a width to the 
mass of the bound-state by the time-energy uncertainty. If this width is a little wider than 
the binding energy A = 2m N — m B , the NN "bound state" can decay into NN (Fig. 6). 
This decay suffers a severe phase space suppression. When the NN bound state is produced 
"on mass shell" {2m^ < M n jj < 2m N + T with T being the width), its decay branching 
fraction to the NN channel is small even if its coupling to NN is strong, i.e., even if g 2 (s) 
in Eq. (fT7|) is large. Consequently, such an NN bound state would appear as a relatively 
narrow meson resonance. On the other hand, when an NN pair is produced above the 
width of the bound state, the bound state can still enhance NN production through the 
small denominator of the resonance propagator. In experiment the transition from "on-shell" 
to "off-shell" occurs continuously above the threshold. The phase-space factor pushes the 
enhancement peak upwards from 2m N to m N jj = 2m N + 0(^J A 2 + T 2 /4). When the phase- 
space factor is removed, the yield curve is expected to behave like g 2 (s)/[(m N jf— ms) 2 +r 2 /4] 
at m N jf > 2m^j. The BES Collaboration [25f first extracted the resonance parameters on 
this assumption when analysis was made only above the pp threshold in J/ip — > ppy. In 
the three-body baryonic B decay, the events of small m N jj receive the SD enhancement, 
as we have argued. The same SD effect would be less prominent in J ftp decay since the 
phase space is much smaller. Consequently the location and the shape of the enhancement 
may not be identical in B and J/ip decays. Magnitude of the net pp enhancement is also 
dependent on dynamical environment of production. Despite such dynamical uncertainties 
we are fairly confident that if enhancement indeed occurs in the region of m N jf > 2m^, 
and if it couples to a state below it, the decay amplitude acquires the negative sign of FSI 
according to the diagram in Fig. 5 and the discussion following Eq. (IT7|) . 



m B 2m N m NN 

FIG. 6: The NN bound state X with mass mg acquires a small width by annihilation decay. 
When the X is produced "on mass shell", only the upper corner (the hatched region) above the 
NN threshold contributes to the decay X — > NN. The broken curve above 2m.N indicates the 
s-wave phase space. The NN production above the X mass shell can still be significant by the 
enhancement due to the X propagator. 

We remark on the coupling between NN and the annihilation channels. In the pp reac- 
tion at the threshold the annihilation cross section is larger than the elastic scattering cross 
section. Can it be compatible with weak coupling between NN and annihilation channels 
? We should first note that the annihilation cross sections fall very rapidly with the inverse 
flux factor l/|v p — v p | above the threshold according to the "l/t>" law of the exothermic 
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reactions. We should also note that the large annihilation cross section is largely due to 
multitude of multi-meson annihilation channels with many different partial waves of sub- 
channels. In contrast, the elastic cross section near the threshold is almost entirely due to 
s-wave scattering. The annihilation cross section may not be so large in many partial-wave 
eigenchannels a little above the threshold. Therefore, the experimentally observed large 
total annihilation cross section is not an outright contradiction with weak coupling between 
NN and meson channels. 

Theorists are not unanimous about existence of the NN bound states and resonances 



19|, |20j. The recent discovery of the state X(1835) in the radiative J/ip decay suggests that 
an NN bound state may exist after all. If X(1835) is indeed an NN bound state in 1 Sq 
or 3 Pq, it is conceivable that a NN bound state exists in the 3 Si, 3 Pi or 1 P\ channel as 
well. Because of negative charge parity, experimental search is harder for 3 Si and 1 P\ in the 
radiative J/ip decay than search of 1 Sq and 3 Pj. Leaving existence of NN bound states as 
an experimental issue still open, we proceed with our hypothesis of the sign flip and study 
the consequences in the pK~ angular correlation. 



IV. NN BOUND STATES AND pK ANGULAR CORRELATION 

The 7T7T7/ resonance X(1835) is the best candidate for the NN bound state. The sharp 
pp threshold enhancement observed in Jip — > ppry first hinted its existence as an iViV bound 



state 251. The mass was deduced at 1859i? (stat)^ 5 (sys) MeV with width < 30 MeV. 
These values are sensitive to the method and assumptions involved in extracting them, e.g. 
rescattering and nonresonant background. They called it X(1859). Two years later the BES 
Collaboration [l3| identified a resonance in the nnr] mass plot and called it X(1835), which 
is presumably the same state as X(1859). It is most likely a state of x So with / = Q[13]. 
Assignment to 3 Po (c and r( in p-wave) of / = cannot be excluded purely experimentally 
though less likely in theory because of the centrifugal repulsion. The width (67.7±20.3±7.7) 
MeV is fairly narrow for its high mass. The upper tail of the width extends beyond the pp 
threshold and contributes to the decay into pp. The BES Collaboration quotes the ratio of 
branching fractions as Br(X(1835) — > pp)/Br(V(1835) — > 7r + 7r~r/) ~ 1/3. In view of the 
tiny pp phase space, we reason that coupling of X(1835) to pp is much stronger than that to 
mesons. For this reason the BES Collaboration suggests that A"(1835) is a likely candidate 
for a molecular or deuteron-like NN bound state. Such a bound state can play a dominant 
role in producing a pp pair in its eigenchannel near the threshold with little annihilation 
into mesons. This is exactly the state that we want for the cause of the sign flip. 

If an NN bound state exists in x So, a bound state may exist in 3 Si as well by the property 
of the meson-exchange force between N and N. If so, the decay amplitudes A K ( 1 So, 0) and 
A K ( 3 Si, 0) flip their signs from the SD ones in Eq. (fT3|) or Eq. (fl4j) to the experimentally 
observed ones in Eq. (flT|) or Eq. ([121) . In this way we would have a chance to obtain 
the observed trend (1 + cos8 p ) 2 for the pK~ angular correlation. Since the B~ — > ppK~ 
amplitudes consist of both 1 = and / = 1 of NN, the sign flip would occur most effectively 
when the 1 = amplitudes dominate over the 1 = 1 amplitudes; 

I^YScO)! » {a^sm |a°( 3 Si,o)| » {A^SiM (20) 

If the 3 Si bound state is in / = 1 instead of / = 0, the second inequality in Eq. (|20!) should 
be reversed in direction. 
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This is our proposal for the resolution of the angular correlation puzzle. As was men- 
tioned below Eq. (THj) . there are other possibilities if a p-wave bound state exists. In those 
cases the s-wave NN bound state should exist only in 1 Sq or 3 Si, not in both. The spin 
splitting is generally weaker in nuclear forces than the orbital-angular-momentum splitting. 
If this prevails in the NN force, the bound states should appear first in the s-wave channels 
and then in the p-wave channels. However, we should keep our mind open to the other 
possibilities of sign reversal in the p-wave amplitudes. In order to make further advance, 
we need to know more about X(1835) and to search for more candidates of the NN bound 
states. In experiment of B meson physics, we are anxious to know the m pp - dependence of the 
angular correlation since it will provide important pieces of information about spin-parity, 
mass and isospin of the bound states. 

As for B°/B — > ppKs, the measurement was made for the oscillating B°-B aver- 



aged over time and therefore no flavor information is available 11]. So long as the pen- 
guin interaction dominates, the time-averaged pKs angular distribution is symmetric under 
cos6*p — > — cos 9 P in general, (cf Eq. (J3]).) Specifically, if we keep only the amplitudes of 
J < 1 and A = in Eq. (Q, the angular distribution is 



dT 



B /B°^ P pK s 



A K CS ,0)\ 2 + \A K ( 3 P ,0)| ; 



+ (|^( 3 5 1 ,0)| 2 + |^( 1 P 1 ,0)| 2 )cos 2 ^]. (21) 



The curve of this angular correlation for pK$ is concave in cos# p . The data[ll| in Fig. 7 
indeed show the tendency of roughly ~ 1 + cos 2 8 P albeit with very large uncertainty. The 
branching fraction was also measured [rD| and its ratio to that of ppK~ is 

Br(B°/B° -> ppK s ) /Bi{B + -> ppK + ) ~ 0.23 (22) 

with roughly ±20% of statistical errors and ±10% of systematic errors. This number would 
be 0.5 if the 1 = amplitudes completely dominates over the 1 = 1 amplitudes in Eq. (!20j) . 
If the 1=1 amplitudes are about 20% of the 1 = amplitudes, however, this ratio 0.23 can 
be reproduced. 



V. Ap CHANNEL 

The threshold enhancement has been observed in other three-body baryonic final states, 
Ap7r + , ppii~ , and AAK + (charge conjugated states combined) as well as in Apj and many 
decay modes of the b —>■ c transition. However, the angular correlation has been measured 
only for a few of them; Ap7r + , Apj, and A c p7r + all in low statistics. Let us look into the 
Ap enhancement observed in ApiT + and Ap7 of B decay and the conjugate, for which the 
dominant process is the penguin decay. 

In the color-dominant process of B — > Ap7r + the spectator <i-quark forms tt + through 
capture by the energetic w-quark that comes directly from the strong penguin interaction. 
Since neither A nor p picks up the slow spectator, their average energies in the B rest frame 
should be comparable. That is, the pir + angular correlation in the Ap rest frame ought to 
be more or less symmetric and flat in cos 9p. This naive SD prediction is in line with 
experiment within large uncertainty: The measured angular correlation [ll[ does not show 
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FIG. 7: The pKg angular correlation in the pp rest frame of B°/B° — ► ppKg [Ref. 10]. 

marked asymmetry nor large variation in cos8p. (See Fig. 8.) The p7r + angular distribution 
can be fitted with any of 1, 1 + |6| 2 sin 2 #p, and (1 — \b\cos9p) 2 with small constant |6|. 
Although the angular distribution is consistent with the SD prediction, one cannot rule 
out a resonance for the Ap enhancement: A resonance of J = leading to dT/dQp ~ 1 is 
certainly acceptable. A resonance of J = 1 is neither ruled out since the flat distribution 
arises with |^ +;0 | 2 + \Al_. \ 2 » |(|^_ ;0 | 2 + \Al +;0 \ 2 ). 
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FIG. 8: The pn + angular correlation in the Ap rest frame of B° — >• Apn + [Ref. 10]. 

Let us turn to Ap-f. In the SD picture the energetic s-quark is emitted against 7 by 
the radiative penguin interaction and becomes the constituent of A. Other quarks and 
antiquarks are produced through strong interaction nearly collinearly against 7 to avoid 
creating a fat gluon. Since the spectator d enters in p, the antiproton is less energetic in the 
rest frame of B° than A. The SD prediction is therefore that p moves along 7 in the Ap rest 
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frame, that is, dT/dQp should rise toward cos 9p = 1, where 9p is the angle between p and 7 
in the Ap rest frame. The observed fry angular correlation clearly shows this trend in line 
with the SD prediction. (See Fig. 9.) 
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FIG. 9: The pj angular correlation in the Ap rest frame of B — > Ajry [Ref. 26]. 



Although the SD prediction is right for A7, a resonance cannot be ruled out for Ap. The 
helicity expansion of the Ap7 amplitude is modified by spin of 7. The helicity of 7 is —1 since 
the s-quark emitted by the penguin interaction is left-handed. Therefore the component of 
total angular momentum along the photon momentum J • p 7 is +1 for Ap in their rest frame. 
The Ap angular correlation is given generally by 



dY 



dm^pdVt % 



= ToE lE<A A; i^(%)r> (A = A F -A A ). (23) 

+Ap7 AaA ? J 



Note here that the first subscript of d{ x (9p) is 1 owing to J z = +1. The Ap enhancement 
cannot be a resonance of J = since J > for d^ x (9p). However, we cannot rule out a 



J = 1 resonance since with 



1 + cos 9p = 2d\ 1 (9 p ) 2 + d\ (9p) 2 , (24) 

the linear correlation ~ 1 + cos6*p arises if the helicity amplitudes for Ap7 happen to obey, 
for instance, 

|^_ ;1 | 2 » 2(|^ +;1 | 2 + |^- ;1 | 2 )' A - + ;i « 0. (25) 

Because of the preferred photon helicity, the P7 angular correlation can be asymmetric under 
cos^p — > — cos#p even with a single partial wave unlike those of pK~ and pn + . Combining 
our observations in B~ — > Ap7r + and B~ — > Ap7 together, we can rule out a J = resonance 
for Ap, but not a J = 1 resonance. However, there is no motivation to call for a Ap resonance 
or bound state at present until we see clear discrepancy with the SD prediction. 
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VI. SUMMARY AND REMARKS 



Three-body baryonic decay modes are favoured over two-body baryonic decay modes 
since a baryon-antibaryon pair may be emitted nearly collinearly in three-body decay. In- 
deed experiment confirms that the invariant mass of the baryon-antibaryon pair is strongly 
enhanced near the threshold in most modes. Although the SD picture appears to describe 
general trends of three-body decays in most cases, we have encountered one clear contradic- 
tion with the SD picture in the angular correlation between p and K ± in — > ppK ± . 

Failure of the SD picture means that some LD effect enters the process of B~ — > ppK~ 
and reverses the angular dependence. We have pointed our finger to the nucleon-antinucleon 
bound states for the cause of sign flip of the SD amplitudes and have given a simple diagram- 
matic explanation for it. To explain the decay angular correlation for B + — > ppK + , we have 
postulated that X(1835) be the 1 S bound state of NN. That is, X(1835) is a molecular 
six-quark state qqqqqq or a deuteron-like state and primarily couples to iViV rather than to 
mesons. Many theorists have made the same or similar proposals on the nature of X(1835), 
with motivations very different from ours. In addition, we need a 3 Si bound state of NN. 
The maximum asymmetry of the angular correlation should occur at the energy where kine- 
matically rising p-wave amplitudes become comparable with the falling s-wave amplitudes. 
The gular dependence will tell us a lot about dynamics near the pp threshold. 

Our argument depends on strong interaction dynamics near the NN threshold that has 
not been proven nor disproven experimentally. Some might feel that we have blown out a 
possible solution to a small puzzle into a farfetched speculation. We cannot counter such 
objections effectively. Our argument presented in this paper is a conjecture or a hypothesis, 
certainly not a theorem. Although our argument is only exploratory and speculative, the 
sign flip by a bound state or a resonance can occur generally and cause failure of the simple 
quark-gluon argument of multi-body B meson decay. 

While our argument of the sign flip is exploratory in nature, we would like to emphasize 
that the partial-wave expansion analysed here will be very useful as a general tool to pene- 
trate into complexity of three-body decay dynamics. To show its usefulness, we have ruled 
out convincingly the resonance hypothesis for the pp threshold enhancement. We have also 
shown that the FSI in s-waves or in p-waves alone should not describe the enhancement 
either. While the angular correlation measurement is not extensive nor accurate enough at 
present, we expect that partial- wave analysis of three-body decay will shed more of new 
light on dynamics of B decay in near future since the data are rapidly accumulating. 
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